Spin-squared Hamiltonian of next-to-leading order gravitational 

interaction 



Jan Steinhoff, Steven Hergt, and Gerhard Schafer 
Theoretisch-Physikalisches Institut, Friedrich- Schiller- Universitat, 
Max-Wien-Platz 1, 07743 Jena, Germany 
(Dated: November 12, 2008) 

Abstract 

The static, i.e., linear momentum independent, part of the next-to-leading order (NLO) gravita- 
tional spin(l)-spin(l) interaction Hamiltonian within the post-Newtonian (PN) approximation is 
calculated from a three-dimensional covariant ansatz for the Hamilton constraint. All coefficients 
in this ansatz can be uniquely fixed for black holes. The resulting Hamiltonian fits into the canon- 
ical formalism of Arnowitt, Deser, and Misner (ADM) and is given in their transverse-traceless 
(ADMTT) gauge. This completes the recent result for the momentum dependent part of the NLO 
spin(l)-spin(l) ADM Hamiltonian for binary black holes (BBH). Thus, all PN NLO effects up to 
quadratic order in spin for BBH are now given in Hamiltonian form in the ADMTT gauge. The 
equations of motion resulting from this Hamiltonian are an important step toward more accurate 
calculations of templates for gravitational waves. 

PACS numbers: 04.25.-g, 04.25.Nx 
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I. INTRODUCTION 



Recently, different methods succeeded in calculating the spin dynamics at higher orders in 
the post-Newtonian (PN) approximation. This is an essential step toward the determination 
of more accurate templates for gravitational waves, to be used in future gravitational wave 
astronomy Among others, some methods were developed to derive Hamiltonians that 
fit into the canonical formalism of Arnowitt, Deser, and Misner (ADM) for nonspinning 
pointlike objects in their transverse-traceless (ADMTT) gauge. The next-to-leading 
order (NLO) spin-orbit (SO) Hamiltonian was calculated in [3fl using the spin equations of 
motion (EOM). The NLO spin(l)-spin(2), or S1S2 in short form, Hamiltonian was obtained 
by using the stress-energy tensor in canonical variables Q|. Various contributions to the 
binary black hole (BBH) spin interaction Hamiltonian, even beyond quadratic order in 
spin, were derived in [5|1 by matching possible static source terms of the field constraints 
to the Kerr metric. In [6j various Hamiltonians, again also beyond quadratic order in spin, 
were calculated by inspecting the (global) Poincare algebra. Astonishingly the momentum 
dependent part of the NLO spin(l)-spin(l), or Sf, Hamiltonian is uniquely fixed up to a 
canonical transformation by imposing Poincare invariance and inspecting possible nonstatic 
source terms for the field constraints in the ADMTT gauge. 

However, the Poincare algebra leaves the static NLO part of the Hamiltonian com- 
pletely unfixed. Though this contribution does not involve many terms, it is difficult to 
calculate the terms that do not contribute to the test-mass limit. This is due to the fact that 
these terms are quadratic in the gravitational coupling. Thus they correspond to a nonlin- 
ear solution of Einsteins field equations for (comparable mass, rotating, and self-interacting) 
BBH within the PN approximation. In our present paper we were able to determine this 
missing piece of the NLO Hamiltonian by a three-dimensional covariant ansatz for the 
source terms of the Hamilton constraint. The arguments that lead to a unique fixation of the 
coefficients are, however, far from being straightforward. Finally, at the end of the paper, a 
comparison with a result from [7] will be made. 

If not otherwise stated, we will use the same notation and conventions as in 4|. 
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II. THE WORLDLINE MODEL 



At higher orders the effective description of inspiraling compact objects crucially depends 
on effects arising from their finite size. Calculations can be greatly simplified if one goes over 
to pointlike objects, or worldlines, from the very beginning. It is thus a good idea to model 
a compact object by a worldline with an appropriate stress-energy tensor and EOM. Dixon 
gave a definition for such a stress-energy tensor as a distribution on the tangent spaces 
of the worldline in terms of four- dimensional covariant multipole moments, see, e.g., Q]. 
These multipole moments represent the finite size of the compact object. The corresponding 
distribution on the manifold was also given, see also as well as the EOM for spin and 
center-of-mass, including all multipole corrections. 

From a strict mathematical point of view a distributional stress-energy tensor makes 
no sense in general relativity. For short, a strict mathematical definition of the product 
of distributions does not exist, but is required due to the nonlinearity of Einstein's field 
equations. This is similar to quantum field theories, like quantum electrodynamics (QED), 
where nonlinearities arise from quantum corrections. However, like in QED, this problem 
can be overcome by applying a regularization and renormalization program when solving 
the field equations. Also notice that, due to the lack of a strict mathematical definition 
for the metric on the worldline, the EOM are stated to be valid only for test bodies by 
Dixon. By using the regularized metric on the worldline, the EOM are now defined for self- 
interacting bodies, too. The state of the art regularization techniques in the ADM formalism 
are reviewed in [lcj] and have already shown to give correct results at least up to quartic 
nonlinear order in the gravitational constant. In this paper we have applied the Riesz kernel 
method. Notice that in the ADMTT gauge no redefinition or renormalization of worldline 
coordinates is needed, at least up to the 3.5PN order for point-mass systems, in contrast to 
the harmonic coordinates approach 

For compact objects at a given PN order, only a finite number of multipole moments 
have to be considered. In order to get the full NLO dynamics one has to add quadrupole 
terms to the pole-dipole expressions for stress-energy tensor and EOM. The stress-energy 
tensor density in terms of delta distributions defined on the manifold must be of the form 



3 



C^T^ = / dr 



t M1/ £(4) + (^ Q (5(4))||q 

The quantities t^ " = t VIJ "" depend only on the 4-velocity given by = where z^(r) 
is the parametrization of the worldline in terms of its proper time r, the spin tensor, and 
quadrupole expressions. Notice that, in general, the quadrupole expressions include not 
only the mass-quadrupole moment, but also the flow-quadrupole moment and the stress- 
quadrupole moment, see, e.g., [l3j|. For the pole-dipole particle t^ uaf3 clearly is zero. In 
contrast to the stress-energy tensor of pole-dipole particles, the Riemann tensor can occur 
at the quadrupole level. Therefore the source terms of the constraints are not assumed to 



be covariant generalizations of flat space expressions in the next section, in contrast to {4]. 
Unfortunately, it is not straightforward to read off explicit expressions for the t^ u "' from 
the definitions in 0, ^3], as they are defined there by linear maps of test functions, also 
involving many implicit definitions, instead of delta distributions. On the other hand, one 
can get explicit expressions for t MiA " by taking (TjQ) as an ansatz and evaluate T'T. = using 
Tulczyjew's theorems [12I. While this is a straightforward calculation, however, it is also 
rather long. 

Further the quadrupole expressions must be fixed. As the quadrupole deformation shall 
completely be generated by the spin here, these expressions are given by a four- dimensional 
covariant ansatz with scalar coefficients in terms of the spin tensor and, in general, also 
the 4-velocity. This ansatz can, of course, be restricted to terms quadratic in spin and to 
terms that contribute to the PN order in question. The coefficients must be fixed by further 
considerations, e.g., matching to the Kerr metric as performed in the next section. 

Even after the four- dimensional stress-energy tensor has been explicitly determined by 
the procedures described above, one still has to follow the approach described in 4| in order 
to get a Hamiltonian in the ADMTT gauge. However, we are only interested in the missing 
static part of the Hamiltonian in the present paper. In the next section we show that the 
three-dimensional covariant source terms in the Hamilton constraint can all be fixed in the 
static case, which is all we need here. 

Important for the ansatz in the next section is the fact that terms like S^ a , where 
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S^ v = —S Ufl is the spin tensor, can only occur as S^u™ in the stress-energy tensor. The 
reason is that the spin tensor is an element of the tangent spaces on the worldline only. Thus 
only the total covariant derivative of the spin tensor with respect to r, i.e., S^^u", is defined 
in this model. However, this total derivative can be eliminated by repeated applying the 
EOM and neglecting higher spin terms. Finally the spin tensor will always be multiplied with 
the delta distribution in the stress-energy tensor directly, i.e., with no derivative operators 
in between. The same applies also to the quadrupole tensor and the 4-velocity or linear 
momentum. 

Notice that in 4| we redefined the center, i.e., we were going over to a new worldline, 
by Lie shifting the source terms of the constraints. For this purpose the linear momentum, 
and in general also spin and quadrupole tensors, had to be defined as a field off the original 
worldline. In 4J we defined this linear momentum field in between the worldlines to be 
parallel transported without rotation, because only with this definition the Poincare algebra 
was fulfilled. Incidentally, this leads to source terms on the new worldline that still fulfill 
the requirements of the last paragraph. On the other hand, one can conclude that one has 
to define the linear momentum between the world lines such that these requirements are 
fulfilled on the new worldline, too. This again leads to the result in [4J without referring to 
Poincare invariance. 

III. ANSATZ FOR THE SOURCE TERMS 

There are still some pathologies that can occur when deriving an ADM Hamiltonian by a 
(3+l)-splitting from the four-dimensional stress-energy tensor, e.g., the source terms of the 
Hamilton constraint 7Y matter = N^/^qT 00 or the momentum constraint 7Y™ attcr = yJ^gT\ 



might depend on lapse or shift, see [4|. Because of the presence of np = —^Jm 2 + Y^PiPj , 
where 7 y denotes the inverse metric to the three-dimensional metric gij, in every source 
term of the Hamilton constraint linear in spin, 7Y matter is the sum of static and nonstatic 
parts only after the PN expansion. This will be true for, at least, most of the quadratic 
source terms, too. However, it was shown in 5| that source terms compatible with the ADM 
formalism sufficient for the nonstatic NLO Hamiltonian do exist, which strongly suggests 
that also in the static case no pathologies occur (or can be overcome, as we only need an 
expression for 7Y matter linear in the metric here, see below). This includes the existence of 



5 



a three-dimensional Euclidean spin vector with constant length as well as the existence 
of a constant masslike parameter m. As Su\ must be given in an Euclidean basis, it can 
be related to a spin tensor in a coordinate basis with the help of a dreibein by 
Sij = ei(k) e j(i) e kimS(m), where tkim is the usual Levi-Civita tensor in flat space, see \&\. The 
dreibein as a function of the metric is just e^u\ = ifj 2 5ij because the metric can be taken as 
conformally flat, gij = ifj A 5ij, in our approximation. 

Our ansatz at quadratic order in spin will be expressed in terms of the three-dimensional 
covariant expressions I[ 3 (mass-quadrupole moment of object 1) and Si given by 

7f = f VS™^™^ + ^Sh ij , (2) 
2S\ = -f^SujSm = const . (3) 

The relation to Qjf and a\ used in |6| is T[ 3 = vn\Q % l and S\ = m\a\. As argued in the last 
section, terms like Suj-f. or Ii 5\-k cannot emerge when coming from the four- dimensional 
worldline model. The most general covariant expression for 7^ matter in the static case then is 

mi ' 3 rrii 



Here Kij and R are the three-dimensional Ricci tensor and scalar, respectively. The q are 
some constants that must be fixed by additional considerations, like matching to the Kerr 
metric. The noncovariant terms are due to the transition from three-dimensional covariant 
linear momentum pi to canonical linear momentum Pj given by 

Pi -Pi- \gin lm i j \Jki + o(p 2 ) + o(s 2 ) (5) 

in Eq. (4.14) of Ref. |4j. Thus the source terms are indeed covariant when the point-mass 
and linear-in-spin terms depending on the (noncovariant) canonical linear momentum are 
added. 

No new static source terms at quadratic order in spin arise in the momentum constraint. 
All we need in the present paper is 

1 



-i/matter 

static 



^ mk S ik 5 + O (S s ) , (6) 
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sec |4. 5]. When adding the point-mass term P,5 one again gets a three-dimensional covariant 
quantity. Though this source term is only linear in spin, it contributes to the quadratic order 
in spin due to the nonlinearity of Einstein's field equations. As already noted in Q| , this is 
the only linear source terms that contributes to the S^-Hamiltonian. 

It is straightforward to calculate the three-dimensional metric ^ arising from this source 
within the PN approximation. Comparing with the three-dimensional Kerr metric in 
ADMTT coordinates from 5] one immediately gets c\ = —\. At a first look it seems 
to be impossible to fixate C2 by comparing to the metric of a single black hole because in 
this case this term does not contribute to the three-dimensional metric. However, notice 
that Tij can be calculated from 7-^ mattcr using 

N ^ S[J d 3 x(NH matter - A^ matter )] , N 

2 ^ Tij = 6^ • (7) 

(Recall that Sf m ) is constant under the variation on the right-hand side, not Sij.) But the 
source term of the lapse function is proportional to NSij^^Tki, and it turns out that 
C2 contributes to the lapse function of a single black hole. Comparing the resulting lapse 
function with the one from we finally get C2 = 0. Incidentally, only this choice for c-i leads 
to the well-known leading order (LO) Sf-Hamiltonian, which is another argument for this 
choice. The last missing coefficients C3 and C4 can still not be fixed, neither by considering 
the three-dimensional metric nor by investigating the lapse function. However, they can be 
set to zero, as they do surprisingly not contribute to the Hamiltonian for BBH at our order 
of approximation. Thus our source terms finally turn out to be just covariant generalizations 
of flat space expressions (using minimal coupling), which was already found to be true for a 
large class of matter couplings to gravity within the ADM formalism in [151 ] . 

Notice that the arguments from the four-dimensional worldline model given in the last 
section led to important and not quite obvious restrictions for the ansatz in this section. 
Without these restrictions a fixation of all coefficients by matching to a single black hole 
would not be possible. Further notice that the source term with coefficient c\ is only needed 
up to linear terms in the metric for the NLO S^-Hamiltonian, because it is a total divergence 
of a vectorial density. This also justifies the validity of the conformal flat approximation for 
the dreibein e^. 
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IV. RESULTING HAMILTONIAN 



The contributions to the Hamiltonian up to NLO quadratic-in-spin terms in the PN 
approximation for BBH can be written as 

H - U i IJLO , jjLO , rjLO 



tt I zjLO I ttLO i rrLO i rrLO 
PM + ^SO + H S 1 S 2 + H S'i + H S 2 2 

) , o-NL 



NLO 



, rrNLO , rr 
+ ^SiS 2 + jH ' 



NLO 



NLO 



ifpM is the point-mass (PM) ADM Hamiltonian known up to 3.5PN, see, e. 
contributions are well-known, see, e.g., [17]. H^q° was recently found in 
4j, also see [18]. 



(8) 

Q. The LO 
and H^o in 



The up to now missing NLO Hamiltonian finally results from the quadratic-in-spin 

r 

static source terms derived in the last section, the nonstatic part of the Hamiltonian from 6] 
and the static linear-in-spin source term of the momentum constraint from [4j. As already 
noted in [4|, the momentum constraint source term linear in spin gives a contribution via 



the integral — ^ / d 3 x < 
is 



(2) (^(3)) 2 to the S^-order (tt 1 ^ is linear in spin). The result for 



rrNLO 



G 



12 



m 2 
4mf 



:ivsi) 2 + 



3m? 



+ 



+ 



3m 2 
Am\ 

3 
Am\ 

3 

2m| 
3 



Pi ■ n l2 ) ( 
Pi • P 2 ) s\ 



5m'{ 

Si ■ n l2 ) 
9 



3m 2 
8mf 



P 2 (5! ■ nia) 



(Pi • Si 



P^ 2 + 



4mi?7T,2 

Px ■ P 2 ) («! • n 12 ) 2 
3 
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P 2 (Si • n 



12 ) 



P x • n 12 ) (P 2 • 50 (5i • n 12 ) + 

15 



m; 



P 2 • n 12 ) (Pi • SO (Si ■ n 12 ) 



4m| 



Pi ■ m 2 ) (P 2 ■ n 12 ) Sf 



4mf 



(Pi ■ n 12 ) (P 2 ■ n 12 ) (Si • n 12 ) 



G 2 m 2 



12 



- ( Si.ni 2 ) 2 --S? + ^-, 



Si n 



12; 



3m 2 
mi 



s 2 



(9) 



Here we have ri 2 = \zi — z 2 \ and rii 2 ri 2 = zi — z 2 , where Zi and z 2 are the three-dimensional 
canonical position vectors of the black holes. results from if<X LO by an exchange of the 

particle labels. This Hamiltonian has the correct test-mass limit, presented, in the ADMTT 
gauge, in Eq. (5.3) of Ref. 

The time evolution of an arbitrary phase space function A, i.e., A = {A, H}, including 
the EOM, follows from the simple, fully reduced equal-time Poisson brackets (with particle 
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labels suppressed) {z\Pj} = 8ij, {*%),<%)} = £ijkS(k)> zero otherwise, where z\ Pj and S^) 
are the components of the vectors z, P, and S, respectively. 

The potential in 7| is not given in reduced spin variables, i.e., with inserted appropriate 
noncovariant spin supplementary condition. The best way to compare is thus by checking 
their EOM for the spin with constant Euclidean length, their Eq. (60). We could not 
find agreement with the EOM for spin resulting from our Hamiltonian H^°. Formally, 



this discrepancy can be attributed to aim in Eq. (61) of Ref. 7]. Through the formal 
replacement 

3 

oi(i) -> Oi(i) + -z(Pi x Si) (10) 

ft v*i 

agreement with our result is achieved. The difference of the precessional frequencies then 
reads 



3 

1 



NW SHS _ ~ (» q \ 

w s 2 _ w s a — 7^3 [Pi x ^ i Pi 



2mf 



(P! X Si) X Pl (11) 



+ A (Pi X 5 l) X Pl 



nxi 



(The index "SHS" refers to our expressions and "NW" to the ones in Ref. Q.) Neglecting 
time derivatives of Si, which obviously are of higher order, this difference can be written as 
a total time derivative and is thus due to a different definition of the spin variables, i.e., 

S NW = S 3HS 3 s ]xSi 

, 2 "" (12) 



2mf 



[(Pi x Si) x p x ] x Si . 



Obviously, the total time derivative of (Tl2|) leads to ffTTl) . Notice that all possible canonical 
transformations can be written as spin redefinitions |3J. According to Eq. (63) or (23) in 18], 
the disagreement given in Eq. (Ill) I) is most likely due to the specific choice of the local frame 
in 0], and thus suggests to consider a redefinition of it. Notice that such a redefinition of the 
local frame also affects Eqs. (24) and (25) in 18J], i.e., it leads to further O (G 2 ) contributions 
to the spin EOM via the algebra. [The suspected redefinition can thus not directly be read 
off from the formal replacement in Eq. (ITU1) .] 

Furthermore, a precession equation derived from a Hamiltonian and standard canoni- 
cal Poisson brackets for the spin variables has to have a special structure. Writing the 
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Hamiltonian as H = ttijS^S^ the spin EOM is given by 

<%) = Eijkifyl + ^lj)S(k)S(i) . (13) 

The precessional frequency reads Uj = (Qji + Qij)S(i). Obviously, Qij is symmetrized in the 
EOM but Eq. (60) of Ref. |7j does not fit into this scheme. However, it can be achieved by 
a spin redefinition which happens implicitly in our Eq. ffl2|) . 

Many of the terms in are proportional to and do therefore not contribute to 

the spin EOM. A comparison of the center-of-mass EOM should thus be envisaged in the 
future. This might also help to clarify the disagreement in the spin EOM. It should be 
noted that our paper goes beyond [7] in the sense that it provides both the spin precession 
and center-of-mass EOM in a condensed form via a Hamiltonian and fully reduced Poisson 
brackets. 
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